Wave-train Induced Unpinning of Weakly Anchored Vortices in Excitable Media 
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A free vortex in excitable media can be displaced and removed by a wave-train. However, simple 
physical arguments suggest that vortices anchored to large inexcitable obstacles cannot be removed 
similarly. We show that unpinning of vortices attached to obstacles smaller than the core radius 
of the free vortex is possible through pacing. The wave-train frequency necessary for unpinning 
increases with the obstacle size and we present a geometric explanation of this dependence. Our 
model-independent results suggest that decreasing excitability of the medium can facilitate pacing- 
induced removal of vortices in cardiac tissue. 
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PACS numbers: 87.19.Hh,05.45.-a,87.19.1p,05.45.Gg 

Rotating spiral waves of propagating excitation char- 
acterize the disruption of ordered behavior in excitable 
media describing a broad class of physical, chemical and 
biological systems [T]. In the heart, spiral waves of elec- 
trical activity have been associated with life-threatening 
arrhythmias [l-Q, i-e., breakdown of the normal rhyth- 
mic pumping action of the heart. Controlling such spa- 
tial patterns with low-amplitude external perturbation is 
a problem of fundamental interest with significant 
implications for the clinical treatment of cardiac arrhyth- 
mias [10]. 

In a homogeneous active medium, a spiral wave can be 
controlled by a wave train, induced by periodic stimula- 
tion from a local source (pacing) [^]. If the frequency 
of stimulation is higher than that of the spiral wave, 
the wave train induces the spiral to drift. In a finite 
medium, the vortex is eventually driven to the boundary 
and thereby eliminated from the system [lll-[l3|. Inho- 
mogeneities in the medium, such as inexcitable obsta- 
cles, can anchor the spiral wave preventing its removal 
by a stimulated wave-train [1^ . This mechanism is anal- 
ogous to pinning of vortices in disordered superconduc- 
tors [15|,ll6|. In the heart, obstacles such as blood vessels 
or scar tissue, can play the role of pinning centres [l7| , 
leading to anatomical reentry, the sustained periodic ex- 
citation of the region around the obstacle. 

In the immediate neighborhood of the obstacle, pinned 
vortices are qualitatively equivalent to waves circulating 
in a one-dimensional ring. They can be removed by ex- 
ternal stimulation provided the electrode is located on 
the reentrant circuit, i.e., the closed path of the vortex 
around the obstacle, and the stimulus is delivered within 
a narrow time interval [l8|. However, for the more gen- 
eral situation of pacing waves generated far away from 
the reentrant circuit, a classical result due to Wiener 



and Rosenblueth (WR) states that, all waves circulat- 
ing around such obstacles are created or annihilated in 
pairs (see Ref. [19], in particular, pp. 216-224). This im- 
plies that it is impossible to unpin the spiral wave by a 
stimulated wave train. 

In this paper, we demonstrate that the WR mecha- 
nism for the failure of pacing in unpinning spiral waves 
is valid only when the radius of the free spiral core (i.e., 
the closed trajectory of the spiral tip defined as a phase 
singularity [20]) is small compared to the size of the ob- 
stacle. We elucidate the transition between the case of a 
free vortex and one attached to a large obstacle by sys- 
tematically decreasing the core radius of the free spiral, 
Rfs^ relative to the obstacle size, Robst^ by increasing the 
excitability of the medium. Our main result is that an 
anchored rotating wave can be removed by a stimulated 
wave train provided Rps > Robst- 

To illustrate our arguments, we use the simple model 
of excitable media introduced in Ref. [21], described by 
an excitatory (u) and a recovery (v) variable: 



dtu = -u{l — 
dtv = {u-v), 



(v + -)] + V^ii, 
a 



(1) 



where, a and b are parameters describing the kinetics. 
The relative time-scale e between the local dynamics of 
u and V is set to 0.02. We discretize the system on a 
square spatial grid of size L x with a lattice spacing 
of Ax = 0.25 and time step of At = 0.01. For our sim- 
ulations L = 200. We solve Eq. [1] using forward Euler 
scheme with a standard nine-point stencil for the Lapla- 
cian. No flux boundary conditions are implemented at 
the edges of the simulation domain. An obstacle is im- 
plemented by introducing a circular region of radius Robst 
in the center of simulation domain inside which diffusion 
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FIG. 1: (a) Wave 5*0, pinned to an obstacle (shaded), ro- 
tates counterclockwise; wave 1 is the first pacing wave, (b) 
Wave 1 hits the obstacle, and separates into a wave rotating 
counterclockwise (la) and a wave rotating clockwise (16). (c) 
Waves 5*0 and \h collide and merge leaving only one rotating 
wave la denoted S\ hereafter, (d) The wave resulting from 
the merging of 5*0 and \h leaves the system. The interaction 
between the following pacing wave, 2 and 5*1, is similar to that 
shown in (a-c). Thus, the pinned vortex persists. Numerical 
simulation of the Barkley model with parameters: a = 0.9, 
h — 0.17; the pacing period is Tp = 6.7 and the radius of the 
obstacle is = 6.5. 



is absent. Pacing is delivered by setting the value of u to 
= 0.9 in a region of 6 x 3 points at the center of the up- 
per boundary of the simulation domain. The maximum 
pacing frequency is limited by the refractory period^ Tref , 
the duration for which stimulation of an excited region 
does not induce a response. 

When the obstacle size is large relative to the core ra- 
dius of the free spiral, Rfs^ the failure of a wave train 
in unpinning the vortex is illustrated in Fig. [TJ Initially, 
the spiral wave 5*0 rotates counterclockwise around the 
obstacle. During the interaction with pacing waves, the 
number of waves attached to the obstacle can change due 
to two possible processes (see Ref. [19], p. 216 and 220). 
First, when the pacing wave reaches the obstacle, it splits 
into two oppositely rotating waves: one clockwise and the 
other counterclockwise. Second, collision between two ro- 
tating waves, as seen in Fig.[T]^c), results in the annihila- 
tion of a pair of counterclockwise and clockwise waves. In 
both cases, the number of waves rotating counterclock- 
wise is always larger than the number rotating clockwise 
by 1. Thus, in addition to conservation of total topolog- 
ical charge (i.e., sum of the individual chiralities, +1 or 
— 1) for all spiral waves in a medium [2o|, [22|, topolog- 
ical charge around the obstacle also appears to be con- 
served. However, in the limiting case of infinitesimally 
small obstacle corresponding to a free vortex, a stimu- 
lated wave-train with frequency higher than that of the 
spiral wave will always succeed in displacing the latter, 
eventually removing it from a finite medium. Thus, there 
is a transition from failure to successful pacing as Robst 
is reduced relative to Rfs- 

The primary fact responsible for this transition is that 
the spiral wave is no longer in physical contact with an 
obstacle of size smaller than Rf s 11 contrary to the 
fundamental assumption of Ref. [l9|. Fig. [2] shows an 
explicit example of successful detachment of a pinned 
wave from the obstacle boundary, where the core radius 
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FIG. 2: Lowering excitability results in successful detachment 
of pinned vortex by pacing. 5*0 is a rotating wave whose core 
(dashed line) is larger than the pinning center (shaded), (a-c) 
are topologically as in Fig. 1. (d) A wavelet Ic is produced 
after collision of waves 5*0 and lb, in contrast with Fig. 1(d). 
(e) The wavelet Ic collides with la and the resulting wave 5*1 
is displaced away from the obstacle, (f) Subsequent pacing 
induces drift of the spiral wave 5*1 to the boundary, eventu- 
ally removing it from the medium. The parameters are as 
in Fig. [U except for a = 0.895 and b = 0.1725, resulting in 
increasing the vortex core size. 



of a free spiral in the medium is made larger than Robst 
by diminishing the excitability of the system. 

The possibility of unpinning the wave in Fig. [2] can 
be traced to the following fact: the collision between 
and the pacing wave-branch lb occurs a small distance 
away from the obstacle boundary and does not result in 
complete annihilation of both waves. A small fragment 
Ic survives in the spatial interval between the collision 
point and the obstacle [Fig.[2fd)]. If the tip of is close 
to the obstacle, the fragment Ic is small, and rapidly 
shrinks and disappears. However, if the gap between the 
reentrant wave tip and the obstacle is large at the colli- 
sion point, such that the size of Ic is larger than a critical 
value In^ the fragment can survive. As Ic propagates fur- 
ther away, it collides with the pacing wave la and forms 
a new broken wave Si that is completely detached from 
the obstacle. Interaction with successive pacing waves 
progressively pushes the vortex further away from the 
obstacle, and eventually from a finite medium. The dif- 
ference between the number of spirals rotating counter- 
clockwise and clockwise around the obstacle changes from 
1 initially (Fig. [21 a), to in Fig. [2fe), contrary to what 
happens for a larger obstacle (Fig. [1]). The absence of 
topological charge conservation for waves rotating around 
a smaller obstacle underHnes the breakdown of the fun- 
damental assumption behind the WR argument for why 
pacing cannot detach pinned waves. The unpinned wave 
is subsequently driven outside the system boundaries by 
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FIG. 3: (a) Parameter space of the Barkley model. Unpin- 
ning is possible in the shaded portion of the SW region, which 
exhibits persistent spiral waves. The thick line indicates the 
boundary with the SE region, where spirals cannot form. The 
domain where unpinning is possible shrinks with increasing 
size of the pinning center, the three dashed lines correspond- 
ing to Rohst — 0, i.e., no obstacle (square), Rohst — 1-25 (plus) 
and Rohst = 6.5 (diamond), (b) Radius Rfs of the free spi- 
ral and the maximum obstacle radius R^^t from which wave 
trains can unpin vortices, as a function of the distance d from 
the SE-SW boundary, along the dot-dashed line indicated in 
(a). Note that Rfs > R^^t^ ^ind both increase with decreas- 
ing d. [In (a), NW (BI) indicates the parameters for which 
steady waves are absent (the medium is bistable).] 



pacing (Fig. [21 f), thus eventually also reducing the total 
topological charge of the finite medium to 0. 

The relative size of the obstacle, compared to the free 
spiral core, is the key parameter that decides whether a 
pinned reentrant wave can be removed or not. Indeed, 
the radius of the free spiral core in the successful case, 
Rfs = 9-05 (Fig. [2]) is significantly larger than in the 
unsuccessful one, Rps = 5.80 (Fig. [T]). It is further con- 
firmed by a detailed numerical study of the interaction 
between a pacing wave train and a pinned spiral over the 
(a, b) parameter space of the Barkley model. As shown 
in Fig. [3l^a), the rotating wave anchored to the obsta- 
cle can be removed by pacing only in the neighborhood 
of the sub-excitable (SE) region (using the terminology 
of Ref. [23]), where Rfs diverges [Fig. M^)]. This is 
explained by noting that in the SE regime, the tangen- 
tial velocity of a broken wavefront is negative, thus caus- 
ing the front to shrink and not form a spiral. As we 
approach the regime where spiral waves are persistent 
(SW), the tangential velocity of the wave break gradu- 
ally increases to zero and becomes positive on crossing 
the SE-SW boundary, so that the broken wave front can 
now evolve into a spiral. As Rps increases with decreas- 
ing tangential velocity of the wave front, the spiral core 
becomes large close to the SE region resulting in success- 
ful pacing-induced termination of pinned reentry. 

We observe that there is a maximum radius of the ob- 
stacle {R^^^f) close to Rfs above which pacing is unsuc- 
cessful in detaching the anchored spiral wave [Fig. [3fb)]. 
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FIG. 4: (a) The maximum pacing period T^^^ at which 
unpinning is possible as a function of the obstacle radius Rohst- 
For the parameters a = 1.1323, b = 0.2459 that we have used, 
the maximum radius of obstacle from which depinning can 
occur is R^}Zt — 4. Tfs is period of a free spiral wave and 
Tref is the refractory period. The dashed line indicates the 
prediction from Eq. [2] (b) The wavelet formation mechanism 
leading to the detachment of the pinned vortex (schematic), 
(c-f) Numerical simulation of the Barkley model. S collides 
with wave 1 at point C at an angle 0. The part 16 of the pacing 
wave merges with moving out of the system. The remaining 
part of the pacing wave collides with the obstacle (shaded) 
separating into la and a small wavelet Ic. When the length / 
of wavelet Ic is larger than the critical nucleation length, Ic 
survives and collides with S. This results in unpinning of S. 



Fig. Ill^a) shows that the pacing period for successful 
unpinning from the obstacle is bounded by the refrac- 
tory period {Tref) and a maximum value T^^^ that is 
independent of Rohst for small obstacles. As we ap- 
proach i^^gf , the upper bound sharply decreases, be- 
coming equal to the refractory time at i^^^f , which indi- 
cates that pacing will be unsuccessful in unpinning waves 
attached to obstacles of radii larger than R^^^ . Thus, 
the results shown in Figs, ^h) and HJ^a) demonstrate 
our earlier assertion that pacing induced removal of an- 
chored waves will be possible only when the obstacle is 
smaller than the core radius of the free spiral wave in the 
medium. 

Our numerical results indicate that the maximum pac- 
ing period necessary for detaching a pinned spiral wave is 
a decreasing function of the obstacle size [Fig.HJ^a)]. This 
can be explained semi-quantitatively by the following ge- 
ometric argument, valid when the size of the obstacle is 
small compared to the core size of the spiral, and sup- 
ported by the simulations shown in Fig. HJ^c-f). The tip 
of the spiral S moves along its circular trajectory, shown 
by the broken line in Fig. Hl^b), and interacts with the 
pacing wave coming from the top, represented by a solid 
line. The part 16 of the pacing wave collides with S at 
the point C characterized by an angle that the spiral 
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tip makes with the symmetry axis (i.e., the hne joining 
the centers of the obstacle and spiral core); the resulting 
wave eventually leaves the system [Fig. H^d)]. The re- 
maining section of the pacing wave splits into two waves, 
la and Ic, propagating along either side of the obstacle. 
The wave tip moves approximately in a straight line from 
C, so that the length of the wave Ic at the symmetry axis 
is / = Rfs{^ + cos6) — 2Robst' When the fragment Ic 
is larger than the nucleation size In^ it expands into a 
wave front that reconnects with wave la. This results in 
a displacement of the wave la away from the obstacle, 
leading to unpinning (as in Fig. [2]). For / < Z^, Ic shrinks 
and eventually disappears, resulting in unsuccessful pac- 
ing. 

Thus, the condition for detachment is / > In- The 
length / is a decreasing function of the angle ^, which in 
turn, is a decreasing function of the pacing period, Tp, 
as explained below. The relation between Tp and can 
be established by estimating the time interval for two 
successive collisions of the spiral with the pacing waves. 
From the point of collision C, the pacing wave reaches the 
obstacle after time Ti = (i^i?^ sin — Robst)/^^ and the 
symmetry axis after time T2 = Ti + {RobstTrs /^Rfs)- 
From the symmetry axis, the new reentrant wave S moves 
by an angle {0 + tt) to arrive at C at time T3 = T2 + 
[Tfs{^ + 7r)/27r], where it collides with the next pacing 
wave. Noting that T3 = Tp allows us to implicitly express 
Tp as a function of ^, and thereby, I. The maximum 
pacing period leading to detachment is obtained when 
/ — ^77,7 as. 

Tp - ^(sm^c-/i?) + ^^ + , (2) 

where, Oc = arccos(2/i^ - 1 + [Iu/Rfs]) and fn = 
Robst/RFS- When Robst > R^Tt = ^fs - (^n/2), T^^^ 
has complex values, indicating that for larger obstacles 
the fragment is too small to survive. The nucleation 
length In can thus be estimated from i^^^f , which al- 
lows us, in turn, to determine the dependence of T^^^ 
as a function of Robst from Eq. [2l Fig. SJ^a) shows this to 
be in fair agreement with our numerical simulations. 

We stress that the arguments used here are model inde- 
pendent, and are based only on the property that waves 
in excitable media annihilate on collision. We verified nu- 
merically that wave-train induced unpinning is observed 
also in a more detailed and realistic description of car- 
diac tissue, the Luo-Rudy model ^] [see the movie on 
line], under conditions of reduced excitability. Meander- 
ing, which occurs in the Barkley model at low a, b values 
(Fig. [3l a), does not affect the physical effect discussed 
here. Note that the proposed unpinning mechanism is 
for the case of an obstacle smaller than the vortex core. 
It is possible under certain circumstances to unpin waves 
from obstacles larger than the core because of other ef- 
fects such as the presence of slow conduction regions [25 1 
and nonlinear wave propagation (alternans) [2 61] . 



Our results thus predict that in cardiac tissue, the re- 
moval of spiral waves pinned to a small obstacle by high- 
frequency wave trains is facilitated by decreasing the ex- 
citability of the medium. This is consistent with previ- 
ous experimental results on cardiac preparations using 
Na-channel blockers [l^ and our prediction could be di- 
rectly tested in a similar experimental setup pj|, [iTj . 

In conclusion, we have shown that for a pinned vortex 
interacting with a pacing wave train, unpinning is pos- 
sible when the size of the obstacle is smaller than that 
of the spiral core. The minimum wave train frequency 
necessary for unpinning in the presence of an inexcitable 
obstacle is higher than that for inducing drift in a free 
vortex towards the boundaries, and it increases with the 
size of the pinning center. Our results suggest that low- 
ering the excitability of the medium makes it easier to 
unpin vortices by pacing. 
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